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We study the semileptonic b ^ c decays of the lowest-lying triply-heavy baryons made from b 
and c quarks in the limit mb,mc 3> Aqcd and close to the zero recoil point. The separate heavy 
quark spin symmetries strongly constrain the matrix elements, leading to single form factors for 
ccb — !> ccc, bbc ccb, and bbb — ^ bbc baryon decays. We also study the effects on these systems 
of using a F-shaped confinement potential, as suggested by lattice QCD results for the interaction 
between three static quarks. 
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I. INTRODUCTION 

Triply heavy baryons are systems of great theoretical interest, since they may serve to better understand the 
interaction among heavy quarks in an environment free of valence light quarks. Besides, being baryonic analogues 
of heavy quarkonium, they might yield sharp tests for QCD. Studying these baryons will be also very useful for 
understanding the three quark static potential. 

With no experimental information available on these systems, previous studies have concentrated on their spectrum. 
To our knowledge the first such study was carried out in 1980 [l| using a QCD-motivated bag model (BM). A mass 
formula was derived by Bjorken in Ref. Q providing predictions for the masses that were larger than those found in 
Rcf. In Ref. [2j the possibility for discovery of the ficcc state was also discussed. More recently there has been other 
phcnomenological mass determinations that include nonrelativistic constituent quark model (NRCQM) calculations 0- 
0], the relativistic three quark model (RTQM) evaluation of Ref. [6], or the Regge approach in Ref. More 
fundamental approaches to the subject include the potential nonrelativistic QCD (pNRQCD) studies of Refs. [sl, Q or 
the QCD sum rule (QCDSR) evaluation of Ref. In Ref. [ul the leading order (LO) pNRQCD result of Ref. 
is usecQ, while a mass calculation that includes next-to- next-to-leading order within the same framework has just 
appeared [T^]. The mass of the triply-heavy baryon ri^bb has been also recently calculated in Lattice QCD (LQCD) 
using 2 + 1 flavours of light sea quarks [l^ . 

Triply-charmed baryon production in the e~^e~ reaction was analyzed in Ref. with the result that the predicted 
production rate was very small. Better perspectives for production are expected at the LHC due to its high luminosity. 
First estimates of the cross section production at LHC were evaluated in Refs. [l5l - [T8l |. A recent evaluation [l^ finds 
that around 10'*-10^ events of triply heavy baryons, with ccc and ccb quark content, can be accumulated for 10 fb of 
integrated luminosity. The authors of this latter work conclude that it is quite likely triply heavy baryons would be 
discovered at LHC. With this in mind, study of their properties beyond spectroscopy seems timely. 

In this work, we will study the lowest lying {J'^ = 1/2+, 3/2+) triply heavy baryons composed of b and c quarkfl 
Heavy quark spin symmetry (HQSS) is of particular interest to study these systems. HQSS is an approximate 
symmetry of QCD in the limit mb,mc ^ Aqcd, and has proved to be an extremely useful tool when dealing with 
heavy hadrons [20l - l22| . This symmetry amounts to the decoupling of the heavy quark spins [l^, [2l| . In Ref. [Ij] it is 
argued that this symmetry cannot be considered as asymptotically valid in heavy-heavy states, since the momentum 
exchange between two heavy quarks might be much larger than Apon , in sharp contrast to the situation for heavy-light 
systems. For mesons with two equal-mass heavy quarks, Ref. |23| argues that the hyperfine splitting Am scales as 
Am ~ niQ Q^(mQ) for sufficiently large heavy-quark mass mg and thus asymptotically increases with ttiq. However, 
~ Off (mg) still approaches zero as mg tends to infinity and the hyperfine splitting becomes negligible compared 
to the total mass. Moreover, the linear behavior is estimated in Ref. (23|1 to take over from the ag{mQ) logarithmic 
fall off for heavy-quark masses in the region of 10 GeV. Hence for systems with two heavy quarks, with masses below 
10 GeV, HQSS should still be valid. That being the case, HQSS should be a useful approximate symmetry to address 
baryons made out of three c and/or b heavy quarks, as we aim to do in this work. 

The study of baryons requires the solution of the three-body problem. In the past we have made extensive use of 
HQSS constraints and have developed a simple variational scheme to find masses and wave functions of single [24 1 
and double I25ll h eavv baryons. We have used the resulting wave functions to study their semileptonic 6 — ^ c ^25-28 
and c — > s,(i [29| decays. The separate heavy quark spin symmetries strongly constrain the matrix elements and, in 
the limit mb,rnc Aqcd and close to the zero-recoil point, they lead to single form factors for all these decays. 

Here, we extend our scheme to study triply heavy baryons. We derive for the first time HQSS relations for their 
semileptonic 6 — ?> c decays from which we can make approximate, but model independent, predictions for some decay 
width ratios. We give absolute values of the semileptonic 6 — >■ c decay widths, as well. We also study the effects in 
these baryons of considering a LQCD inspired three-body confinement potential (denoted as F in [i]) instead of the 
commonly used one, obtained from the sum of two-body quark-quark terms. 

II. SPIN SYMMETRY 

We will consider decays induced by the semileptonic weak decay of a 6 quark to a c quark. Near the zero recoil 
point, the velocities of the initial and final baryons are approximately the same. If the momenta of the initial and 
final baryons arc = "mv^ and p'^ = m'v'^ = m'v^ -\- respectively, then k will be a small residual momentum near 
the zero-recoil point. For the initial baryon at rest we have that k ■ v = E' — m! . For a small final momentum this 



^ This coincides with the 1/r, or Coulomb, interaction that comes from one-gluon exchange 

^ In what follows we will denote by S the baryons with spin 1/2, while we will use S* and O* for the spin 3/2 ones. 
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is approximately given by p'"^ /2m' and then is 0(1/™') close to zero recoil. We will work near zero-recoil and thus 
neglect v ■ k below. 

The consequences of spin symmetry for weak matrix elements can be derived using the "trace formalism" [i^, [sO] ■ 
The scheme advocated here is similar to that employed in Ref. [sT] to study the semileptonic be to cc baryon decays. 
To represent baryons with three heavy quarks we will use wave functions comprising tensor products of Dirac matrices 
and spinors. For Q1Q1Q2 baryons containing two heavy quarks Qi and a distinct heavy quark Q2, we have: 



^QiQiQ2 



"Q1Q1Q2 



(1 



a/3 



1 



(1 



-7/x- 



(1 



(1) 
(2) 



al3 



where we have indicated Dirac quark indices a, /3 and 7 explicitly on the right-hand sides and r is a helicity label for 
the baryor0. For the 5* states, ul^{v,r) is a Rarita-Schwinger spinor. For the baryon containing three heavy quarks 
of the same flavour, we use: 



1 



(1 + y) (l-i) 



<(w,r) 



(4) 



These wave functions can be considered as matrix elements of the form {0\Qi^QlpQ2^\BQ-^Q^Q2) where Q'^ — Q^C 
with C the charge-conjugation matrix. In each case we couple two quarks of the same flavour in a symmetric spin-1 
state in the first factor and combine with a spinor for the third quark. Under a Lorentz transformation. A, and heavy 
quark spin transformations 5q, a wave function of the form TapU^y, with 14 — -^{'^'^ + 7^)75W or u^, transforms as: 



TU ^ S{A)rS-\A) S{A)U, TU^Sq.TSI Sq^U. 



(5) 



The Q1Q1Q2 states have normalization Z^WTr(rr), while for QQQ the normalization is UU Tr(rr) + 2U TTU (which 
can be understood by counting quark contractions). We define F = j'^r'^j^ as usual and our spinors satisfy uu — 2m, 
m^'m^ — ~2m where m is the mass of the state. 

We construct amplitudes for semileptonic decays determined by matrix elements of the weak current — C7^(l — 
75)6. The operator cJ^h, where J'^ = 7^(1 — 75), would be invariant under heavy-quark spin transformations if 
transformed as — >■ ScJ^Sl- Thus, we can build matrix elements respecting the heavy quark spin symmetry by 
constructing quantities which would be invariant under the same assumption. We observe that j'' can be rewritten 
as ji" = -b'^^i'il + 75)c'= and note that b^J^f^c", where J"'' = --/''{l + 75), would be invariant if J'^^ ^ SbJ^^'Sl 

For the transitions S^^j, — >■ il*^^, the matrix element respecting heavy quark symmetry is, up to a scalar function of 
the product of velocities, w = v ■ v', 



{ni^,, V, k, r'l j''(0)|s(*) , V, t) = U'{v, fc, r')jm{v, r)Tr[r,,f,r,,J + U\v, fc, r')T,,bTcccJ''U{v, r) 



(6) 



where r and r' are the helicities of the initial and final states, and we use the standard relativistic normalization 
for hadronic states. Terms with a factor of 1/ can be omitted because of the equations of motion (^u = u, = F, 



7/.W' 



0, VfiU^ = 0), while terms with ^ will always lead to contributions proportional to u • fc which is set to at 



the order we are working. We also make use of the exact relation u'^u = and the approximate ones u'^jU — u'v^u, 
u'^^u = 0, and u'l/^^^^u — —u'l/v^^i^u valid close to zero recoil. 



For the transitions .^[^^ — )■ .^[.*^ the matrix element is 



(5^*i^ V, fc, r'|j''(0)|S^;i, V, r) = U'{v, fc, r')n,,J-^T,,^{v, r). 



:(*) 



(7) 



while for the transitions 51 



bbb 



^bbc^ the matrix element respecting heavy quark symmetry now reads 



(S^;^, V, fc, r'\j^{0)\ni,„ V, r) = U'{v, fc, r')J^'U{v, r)TT[T bbbT bbc\ + U'{v, fc, r')TbbbTbbcJ''U{v, r) 



(8) 



Note that the two identical heavy quarks Qi can only be in a symmetric spin 1 state. The structure 

(1 

■7m- 



(3) 

guarantees that the spin of the first two heavy quarks is coupled to 1 (see for instance Refs. [2ll l22l| ). On the other hand, the spin-1/2 



spinor |^-i=(i;'' -l- 7'')75u(d, r)j is discussed in Ref. [32 
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Close to zero recoil, and within the approximations mentioned above, our results for the transition matrix elements, 
apart from irrelevant global phases, are: 

Ecct^^*ccc '^Vu'^u (9) 

Kcb Kcc U'^ni - 15)UX (10) 

5 

Ebbc Eccb -Xu'iY' - g7^75)M (11) 

Ehbc ^ Kcb ^-^^^ 

^Ibc ^ ^ccb --^xw'u'^ (13) 

^Ibc ^ ^*ccb -2xu'^ril - l5)ux (14) 

ni,,^Ebbc 2^u'u^ (15) 

ni,, ^ S*,, - u"ni - 75)^A (16) 

where the factors ri{w), x{^) ^^'^ C(^) the Isgur-Wise functions that depend on w = v ■ v' and that we expect to 
be close to 1 at zero recoil (w = 1). In fact in the limit rric = rrib they would be exactly 1 at zero recoil. To check 
that assertion let us consider an SU{2) symmetry under which the c and b quarks transform as a doublet and the four 
states n*^^, 5^^^, and fllf^^ form a quadruplet. We will consider all heavy quark spins aligned, that is to say we 
will place the four baryons in the state with maximum third component of spin, Jz = +3/2. The fi — component of 
the vector part of the transition operator j^^ would then be = c^6, which is the raising operator in the Fock space 
for this flavour SU{2) symmetry. Assuming this symmetry, and taking into account the state normalization, we will 
have at zero recoil 

V3 = ^^{n:jc^E:,,) = ^/3 77(l) (17) 

2-^(s:,b|ct6|S,V)=2x(l) (18) 

V3=^{El,,\c^b\ni,,)=V3ai) (19) 

In the above equations, the left-most results follow from the c O & SU{2) symmetry, while the right-most ones are 
obtained from Eqs. (HU]), (HI) and From Eqs. we deduce ?7(1) = x(l) = ^(1) = 1- For the actual 

quark masses one expects deviations from this result as a consequence of a mismatch between the initial and final 
wave functions. 



III. DECAY WIDTH FOR A SEMILEPTONIC b^c TRANSITION AND HQSS CONSTRAINTS 

The total decay width for semileptonic b ^ c baryon transitions, is given by 
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r = \Vcb\^^ — / Vw^-ic^^iq)napiv, k) dw (20) 

OTT^ m J 

where \Vcb\ is the modulus of the corresponding Cabibbo-Kobayashi-Maskawa (CKM) matrix element for a 6 — > c 
quark transition, for which we shah use \Vbc\ = 0.0410 33]. Gf = 1.16637(1) x 10"" MeV"^ [33] is the Fermi decay 
constant and q — p~p' . w and q are related hy w = ™ 2mm~'^ • ^he decay, w ranges from w — 1, corresponding to 
zero recoil of the final baryon, to a maximum value given, neglecting the neutrino mass, hy w — Wmax — *" "2mm' ' 
where mi is the final charged lepton mass. Finally £°'^{q) is the leptonic tensor after integrating over the lepton 
momenta and 'Hq^(w, k) is the hadronic tensor. 
The leptonic tensor is given by 

£"'5((Z) = A(g2)g"/J + B(q2)9V (21) 

Q 
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where 



with 



Hi') 



(22) 



The hadronic tensor reads 



I{q') = ^,{q'-mf) 



^±-^Y.{B',v,k,r' \riO)\B,v,r) {B\v,k,r' |/(0)|i?, z;, r )* 



(23) 



(24) 



with J the initial baryon spin. Baryonic states are normahzed such that 

(B, v', r' \B,v,r) ^2E {27rf Srr' S{p- p') 
with E the baryon energy for three-momentum p. 



(25) 



A. HQSS constraints on semileptonic decay widths. 

For large quark masses and near zero recoil we can use the HQSS results in Eqs. ((5])- (fTC)) to approximate the 
product C^f^ Hafi by 

• ccb — > ccc transitions 



1 « 

£"^Ha^ « —r,^mm'{l + w) 



- 3A{q^) + S(g2) 



CCD CCC 



1 



8A{q')w{l + 2w'') + B{q') { - w{l2 + 8w^) + 2^^^^^^§^(20 + 8w^) 



• bbc — > ccb transitions 

^bbc ^ ^ccb 



C'^Hap « -x'mm'<{ - A{q''){34w + 32) + ^(g^) 



(26) 



(27) 



(28) 



^bbc ~^ ^ccb 



3Aiq^) + B{q')[^-^ 



- 1 



(29) 



^bbc ^ ^ccb 



^bbc ^ccb 



iA{q^) + B{q^) 



q^ ) 



- 8A{q^)w{l + 2w'^) + B{q^) ( - w{12 + 8w^) + 2^^^-^^^^^-^ (20 + 8w^) 



(30) 



(31) 



• bbb — > bbc transitions 

^ ^bbb ~^ ^bbc 



3Aiq^) + B{q^)[^-^ 



- 1 



(32) 
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— O* —i -* 
"bbb ^ ^bbc 



r 



(33) 



In the strict near zero recoil approximation, w « 1 or equivalently very close to its maximum value g^axi '^^n. 
approximate 

q2 ql ql 

In addition, A{q^) w -B{q^) near f/j^jax- To the extent that the former approximations are good and further using 

ruB,^, « rriBi^^ ; ms„, « ms-^^ (35) 

we can make approximate, but model independent, predictions for ratios of semileptonic 6 — )■ c decay widths based in 
the above HQSS relations for Cl^Uap. We find 



r(Sfcbc — > Sccb) 

'^^i^bbc '^ccb) 

r(Sb6c ^ s*^f,) 

T('='bbe ~^ "ccfe) 
b^^^bbb ~^ ^bbc) 



1 (36) 

1 (37) 

1 (38) 

1 (39) 

1 (40) 



These relations are similar to the ones we obtained in our former study of doubly heavy baryons 27] and from the 
findings of this latter work, we expect them to hold at the level of 20%. To estimate the decay widths themselves, 
we need to know the Isgur-Wise functions 'q{'w),x{w) and ^{w). In the next section we will use a nonrelativistic 
constituent quark model for this purpose. 



IV. NONRELATIVISTIC QUARK MODEL EVALUATION OF THE ISGUR-WISE FUNCTIONS AND 

DECAY WIDTHS 



In this section we shall obtain, within the nonrelativistic quark model and using the ALl interquark potential of 
Refs. [3, 1131, the wave functions of the heavy baryons involved in this study. With those wave functions we can 
evaluate the Isgur-Wise functions and estimate the baryon semileptonic 6 — >■ c decay widths. 

The wave functions have the general form 

^aia^aa = '^/i/i'^/^/i'^/:,/!' ^^^*(?'l , , ^12) (1/2, 1/2, 1; Si, S2, Si + S2)(l, 1/2, J; Si + S2, S3, M) (41) 

where aj represents the spin (s) , flavour (f ) and colour (c) quantum numbers of the j-th quark. The two flrst quarks 
have the same flavour /i, while the third quark has flavour /i', which could be also the same as the one of the first two. 
eciC2C3/'\/3! is the fully antisymmetric colour wave function and the (ji, j2, j; mi, m2, m) are SU{2) Clebsch-Gordan 
coefficients. J is the total spin of the baryon. As we are interested only in spin 1/2 or 3/2 ground state baryons, the 
total orbital angular momentum is i = 0. Thus, the orbital part of the wave function can only depend on the modulus 
of the relative distances between the quarks. Here we use ri, r2 which are the relative distances between quark three 
and quarks one and two respectively, and ri2 which is the relative distance between the first two quarks. Following our 
works on single and double heavy baryons [l^, we shall use a variational ansatz to solve the three-body problem. 
We write the orbital wave functions as the product of three functions, each one depending on just one of the three 
variables ri,r2,ri2, i.e. 



$(ri,r2,ri2) = (I)hh'{ri)(f>hh'{r2)4'hh{ri2) 



(42) 
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[13] 


[1] 


f2] 


f4] 


m 


[5j 


[6] 




fio] 




Variational Faddeev 


LQCD 


BM 




NRCQM LO pNRQCD NRCQM RTQM 


Regge 


QCDSR 




14398 


14398 


14371 ± 12 


14300 


14760 ± 180 




14370 ± 80 


14834 


14569 


- 


13280 ± 100 




11245 




- 


11200 


11480 ± 120 




11190 ±80 


11554 


11287 


- 


10540 ± 110 




11214 


\V1VJ 










11190 ±80 


11535 


11280 




10300 ± 100 




8046 






8030 


8200 ± 90 




7980 ± 70 


8265 


8025 




7450 ± 160 




8018 


8019 










7980 ± 70 


8245 


8018 




7410 ± 130 


rrr 


4799 


4799 




4790 


4925 ± 90 


4632 


4760 ± 60 


4965 


4803 


4819 ± 7 


4670 ± 150 



TABLE I. Masses (in MeV) of the triply heavy baryons obtained with the ALl potential of Refs. [3, using our variational 
approach. For comparison we also show the results from the Faddeev calculation performed in Ref [3] using the same potential. 
Predicted masses within other theoretical approaches are also compiled. Hyperfine splitting is neglected in [Tll |. 



For each of the functions above we take an expression of the for nfl 

4 

c/,(r) = ^a,e-''?('-+''^)' (ai = 1). (43) 

The variational parameters are fixed by minimizing the energy and the overah normahzation is fixed at the end of the 
calculation. The results we get for the masses are given in Table HI where we also compare them to the ones obtained 
in Ref. j3] using the same potential but solving Faddeev equations. The agreement between the two approaches is 
very good. 

There is a recent estimate [111 of the mass of the triply-heavy baryon hhh obtained in lattice QCD with 2±1 flavours 
of light sea quarks. Our result compares with it rather well. Our predictions are also in a reasonable agreement 
with those obtained within the BM, RTQM and LO pNRQCD evaluations of Refs. @ and [ll|, respectively. 
The QCDSR masses calculated in [10| come out systematically much smaller than ours, while those obtained in the 
NRCQM of Ref. ^ are significantly larger than our predictions. 

Before computing the Isgur-Wise functions that govern the semileptonic decays of the triply baryons, we would 
like to devote a few words to discussing the confinement potential in these systems. In phenomenological constituent 
quark models, such as the ALl potential used here, the confinement potential for baryons is usually obtained from 
the two bod y fo rces that describe the dynamics of each quark pair. However, the lattice QCD simulations carried 
out in Refs. [SSl, |36| seem to indicate that in a three static quark system, confinement is a genuine three-body effect. 
Changes in the masses due to the use of one or other of these approaches are studied in the next subsection. 



A. A-shaped versus y-shaped potential 

LQCD results for the static quark-antiquark ground-state potential [131 are well described by a dependence 

- - + CTr + C (44) 
r 

which contains the sum of the short distance Coulomb one-gluon exchange (OGE) term plus the confining long 
distance flux-tube contribution. Most phenomenological models assume such a dependence and fit the yl, cr and C 
parameters to the meson spectrum. This is for instance the case of the ALl potential that we use. When going to 
the quark-quark sector, a factor of 1/2, assumed to come from an overall colour A • A dependence!!, is added to the 
interaction. The resulting potential in baryons is thus obtained as the sum of two-body terms. For the confining part 
one is summing over the three sides of a triangle with the quarks located at its vertices \<j(r\ +r2 +ri2)/2], leading to 
the name of A-shaped potential (see Figure [1]). This picture works very well from a phenomenological point of view 
and one gets a good description of the light and single heavy baryon spectrum once the parameters have been fixed 
in the corresponding meson sector. 



^ We use four gaussians in the present approach. We have checked that by increasing the number of gaussians, the variational baryon 

masses change in less than 5 MeV. 
^ QCD predicts exactly this colour factor for the OGE term. 
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> / 

FIG. 1. Illustration of the A- and F-shaped confinement potentials. 

As mentioned above, the 3 quark static potential has been directly measured on the lattice in Refs. [S^H^. A good 
fit to the lattice data was obtained assuming a picture, similar to the one described above, in which the potential has 
a short-distance Coulomb OGE part plus a long-distance flux-tube part 

- A3g(— + — + —) +a3qL^in + C3q (45) 

where Lmin is the minimal value of the total length of the colour flux tubes linking the three quarks. The flux tubes 
adopt a Y shape (See Fig.[T]), hence the name F-shaped potential. This is in agreement with the picture that emerges 
from the QCD BM calculations carried out in Ref. Indeed the A and Y nomenclature was already used in this 
pioneering work of 1980. In terms of the ri,r2 and ri2 interquark distances one has 

iiin = \{rl+ rl + rl,) + ^ sj - Xirl rj, rj,) (46) 

when none of the angles of the three quark triangle exceeds 2tt/3 and where A(a, b, c) — a'^ + b'^ + — 2ab — 2ac — 2bc. 
If one of the triangle's angles exceeds 2tt/3 then Lmin is just given by 

imin = + 7-2 + ri2 - max(ri , r2 , ri2 ) (47) 

Comparing this fit with the one for the quark- antiquark potential they found that a^g ~ a, A^g w and Csg w |C. 
Thus, leaving out the confinement piece, one could approximate the 3 quark potential by the sum of three two-body 
quark-quark terms. For the confining part, the sum of three two-body quark-quark terms is always smaller than the 
three body force obtained from lattice QCD data. Actually, one has (ri + r2 + ri2)/2 < Lmin < III + ^2 + ri2)/V^, 
which might induce changes of around 15% at most in this part of the potential. Indeed, in Ref. [35| the lattice data 
were also fitted to the sum of the three quark-quark potential^ and it was found that a slightly larger value for the 
confinement coefficient (0.53 a vs cr/2) was required. The ad-hoc factor of 1/2 introduced in quark potentials when 
going from the mesons to baryons is understood here as a geometrical effect rather than as a colour factor as it is 
usually presented. 

To be more quantitative, we have computed the triply heavy baryon masses also with a K-shaped confinement 
potential. To that end, we have taken the ALl potential used before, and have replaced the a{ri + r2 +7'i2)/2 term by 
cimin- Results are presented in Table |TT1 There we also compare with the masses obtained previously using the ALl 
potential. We see a small increase in the masses of roughly 26, 34, 40 and 48 MeV for the bbb, bbc, ccb and ccc systems 



The fit is worse than that obtained when the functional form of Eq. 1145 l l is used. 
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This work 
A-shaped potential K-shaped potential 



mnj^^ 14398 14424 

m=. 11245 11281 

"fafac 

ms,^^ 11214 11247 

ms;^^ 8046 8087 

ms.^i, 8018 8058 

mnj^^ 4799 4847 



TABLE II. First column: Masses (in MeV) of the triply heavy baryons obtained with the ALl potential of Refs. [1, HJ] using 
our variational approach. Second column: The same by substituting a{ri + r2 + ri2)/2 by aLmin in the ALl potential. 

respectively. Effects here are similar to those due to the hyperfine splitting. Future and precise measurements of the 
masses might help to shed light on the exact nature of the confinement potential in the baryon sector. 

The corrections to the ratios of decay widths, that would be computed in the subsections below, are even smaller, 
as expected from perturbation theory, since changes in the wave functions arise at second order. 

B. Isgur-Wise functions 

To evaluate the Isgur-Wise functions we follow our work in Rcf. 26] and wr 

VE^m = J dVidV2e'^--("=^-^+'"=^~^)/^-[vI-f,.^^(ri,r2,ri2)]* vf^c., (n, r2, na) (48) 

J ^cch ^bbc 

eH(.)(|g|) = J dVid3r2e''-("^^^^+™^^^^)/*^''-[*^,.,(ri,r2,ri2)]* *o-,Jri, ra, m) (50) 

where Mhhh' = f^h + "m-h + ■ In Fig. [2l we display the eight overlap functions obtained from each of the decays 
examined here. We see that as predicted by HQSS in Eqs. they reduce to only three independent ones in 

very good approximation. In the equal mass case they would be equal to one at zero recoil {\q \ = 0). For finite 
masses we see they deviate slightly from 1 at zero recoil owing to the mismatch between the initial and final wave 
functions. Note that w — 1 -\- v ■ k/m! , and thus as w departs from 1, v ■ k/m! increases. To obtain the relations of 
Eq. (|9|)-(fT6|) all V ■ k/m! corrections were neglected. Thus, the Isgur-Wise (overlap) functions depicted in Fig [2] would 
provide a poorer description of the weak transition matrix elements as w deviates from the zero recoil point. The 
largest corrections are expected for the hcc — ccc transitions, related to the ry— type Isgur-Wise functions in Fig [21 
for which Wmax ~ 1.125. For this case at the = end of the phase space, v ■ k/m' becomes of order 1/8. Thus, in 
this region, approximating the full amplitude (weak matrix element) by means of the rj(w) Isgur-Wise function could 
be subject to uncertainties of order of 15-25%. For hhc — >■ bcc or bbb — > hhc transitions, the %— and ^— Isgur-Wise 
functions should provide more accurate estimates of the transition weak matrix elements for the whole available phase 
space, since in those cases Wmax is only about 1.06 and 1.03, respectively. 

In the next subsection, we will make use of the Isgur-Wise functions of Fig [2] to estimate the decay widths. This 
should be quite accurate, even for the bcc — ccc transitions, since, as we will see, the differential decay width 
distribution peaks very close to the zero recoil point and hence far from the end point of the spectrum w — lUmax- 
Indeed, for bcc — > ccc transitions, the distribution takes its maximum value well below w — 1.05 (see left upper panel 
of Fig.El). 

C. Decay widths 

We now use the HQSS approximate expressions in Eqs. (l26)) - (p3)) to estimate the decay widths. The results for the 
differential distributions are shown in Fig. |3] and the integrated decay widths are compiled in Table IIIIl We see the 



Note when the initial baryon is at rest, w = ^ is just a function of \q\. 
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1.025 1.05 1.075 1.1 1.125 1.15 



w 



FIG. 2. Overlap functions for b ^ c semileptonic decays of triply-heavy baryons obtained in a nonrelativistic quark model. 
The functions fall into three families, consistent with heavy quark spin symmetry. 



ratios in Eqs. (|37|) and (|38l) are satisfied at the level of 5.5% and 3.4% respectively whereas the ratios in Eqs. (|36|) . 

and (1^0)) are good only at the level of 20-30%. It is clear the relations in Eqs. can only be approximate. 

First, the strict zero recoil point is forbidden by phase space, and second changes rapidly from its maximum value 
of (m — m')^ at w = 1 to its minimum value of mf at Wmax which makes the approximation in Eq. (j34p not good 
enougtd. 

What one sees when looking at the differential decay widths in Fig.[3]is that these distributions peak in each case in 
the lower part of the allowed w region, about 1.005, 1.009 and 1.025 for bbb, bbc and ccb decays respectively, quite close 
to the zero recoil point. In these circumstances one can relax the strict approximation in Eq. p4p and use instead 

{vqf ^ {vq){v' -q) ^ {v' ■ qf 
qz 

which should be reasonable near the maximum of the differential decay width, since we can still use w « 1. We can 
also use B{q^) k, —A{q^) and the approximate equality of masses in Eq. ([55)) . One can now only make the following 
two model independent predictions 



AT{Ebbc ^ 5ccb) - lOTiEbbc ^ 
which we see are good at the level of 3.4% and 0.25% respectively. 



1 (52) 
1 (53) 



Note, as pointed out in Ref. [27} . the quantities {v ■ q)^/q^, {v' ■ q){v ■ q)/q'^ and {v' ■ q)^ jq^ which are all equal to 1 near zero recoil, 
quickly deviate from 1 because of the (p' factor in the denominator. 



1 





ni , \ , \ , \ 1 

1 1.01 1.02 1.03 

W 

FIG. 3. Estimated ^ differential decay widths in ps~^ for the different transitions considered. 
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7.74 X 10" 
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3.95 X 10" 


^bbb " 




6.34 X 10" 



TABLE III. Estimated decay widths in units of ps ^ . We use |Vi,c| = 0.0410. Similar results are obtained for /xi^^ leptons in 
the final state. 



V. CONCLUSIONS 



We have studied the b ^ c semileptonic decays of the lowest lying triply heavy {Q1Q2Q3, with Qi = b,c) baryons 
in the limit mii,mc ^ Aqcd and close to the zero-recoil point. The separate heavy quark spin symmetries strongly 
constrain the matrix elements, leading to single form factors for all these decays. We have obtained these HQSS 
relations for the first time. Lattice QCD simulations work best near the zero-recoil point and thus are well suited to 
check the validity of our results. 

We have used a NRCQM, adjusted to the meson spectrum, to predict the masses of these triply heavy baryons by 
using a simple variational scheme. Results for masses compare rather well with some previous Faddeev and LQCD 
estimates. We have also obtained masses by using a lattice QCD inspired three-body confinement potential. The 
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variational wave functions have been employed to compute the overlap integrals needed to evaluate the relevant Isgur- 
Wise functions that describe these decays. We have checked that our calculations are consistent with HQSS and have 
used them to estimate the semileptonic decay widths. 

We have in addition made approximate, but model independent, predictions for ratios of semileptonic 6 — >■ c decay 
widths based on the HQSS relations derived here, which we expect to be accurately fulfilled. 
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